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Abstract 

We study the static as well as the glassy or dynamical transition in the mean- 
field p-state Potts glass. By numerical solution of the saddle point equations we 
investigate the static and the dynamical transition for all values of p in the non- 
perturbative regime p > 4. The static and dynamical Edwards- Anderson parameter 
increase with p logarithmically. This makes the glassy transition temperature lie 
very close to the static one. We compare the main predictions of the theory with the 
numerical simulations. 
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1 Introduction 



This work is devoted to the study of the glassy properties of the mean-field Potts glass. Very 
recently there have been new developments in the spin-glass theory concerning frustrated 
mean-field models without explicit disorder [1, 2, 3, 4, 5, 6, 7]. It has been shown that 
these systems do have a glassy transition temperature below which thermal fluctuations 
are very small and dynamical relaxations are very slow. Even though these results are not 
new in the context of disordered systems it is much interesting to know that non disordered 
models also share these properties. 

The purpose of this work is to study the glassy behavior of a disordered spin glass. 
In general, these systems have a static transition Trsb where replica symmetry is broken. 
The breaking of the replica symmetry can occur in two ways. There can be a continuous 
breaking pattern (as happens in the case of the Sherrington-Kirkpatrick (SK) model [8]) or 
there can be a one step breaking of the replica symmetry (as happens in p-spin models with 
p > 2 [9]). Also one can find intermediate phases where there is a pattern with one step of 
breaking superimposed to a region with continuous breaking (as happens in p-spin models 
or Potts models at low enough temperatures). The breaking pattern is fully described by 
the order parameter q(x) which is a function defined in the interval (0, 1) [10]. 

All these transitions are continuous from the thermodynamical point of view, i.e. there 
is no latent heat. But in some cases they can be first-order in the order parameter. This can 
occur because the thermodynamic potentials are continuous (they are usually expressed as 
integrals of the order parameter function q(x)). 

The purpose of this work is the study of a disordered model with a discontinuous tran- 
sition in the order parameter. These systems generally have a temperature Tq where a 
dynamic instability appears. This temperature is called the glass temperature and is higher 
than the transition Trsb where the replica symmetry breaks. The first observation of this 
type was due to Kirkpatrick and Thirumalai who solved the off-equilibrium dynamics for 
the p-spin model above the glass temperature [11]. Subsequently, Kirkpatrick, Thirumalai 
and Wolynes studied the Potts mean-field glass reaching similar conclusions [12, 13]. Sim- 
ilar results were obtained in case of the p-spin spherical spin glass by Crisanti, Horner and 
Sommers [14]. Below the glass transition it has been shown by Cugliandolo and Kurchan 
[15] that the energy of the p-spin spherical spin glass model in the low temperature phase 
is higher than that predicted by the statics. For times larger than a time scale (which 
diverges exponentially with the size of the system) it is expected that the energy of the 
system will relax to its equilibrium value. How fast this time scale grows with the size of 
the system depends on particular features of the glass transition like the discontinuity in 
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the Ed wards- Anderson parameter qQ. 

In order to investigate the glassy behavior of a disordered model we have decided to 
study the infinite-ranged Potts glass model. The reason is threefold. On the one hand, 
in the Potts model p is a tunning parameter for the magnitude of the static and the 
dynamical transition. On the other hand, the Potts model is amenable of numerical tests 
while other models like the p-spin model (Ising or spherical) are time consuming which 
makes numerical simulations practically impossible for p > 3. The situation is different in 
case of the random orthogonal model [2] where the replica theory predicts the existence of 
a glassy phase in good agreement with the numerical simulations. Finally, the Potts glass 
model lacks the reflection symmetry <7j — > — <j{ of some other models. This makes it more 
similar to real structural glasses. 

As we will see later the Potts glass model present one (not serious) drawback. This is 
that for p > 2 the systems order ferromagnetically at low enough temperatures. In order 
to investigate the spin-glass behavior it is necessary to introduce an additional antiferro- 
magnetic coupling constant. It is in these conditions that we have investigated the glassy 
features of the discontinuous transition. The research of a glass transition in presence of 
ferromagnetic order remains an interesting open problem. 

This work completely solves the replica equations for the Potts model for arbitrary 
number of states p. We will be able to exactly compute the static and the dynamical 
transition and we will compare the predictions with the numerical simulations. 

We will see that a complete dynamical freezing never takes place, even for very large 
values of p. For generic p there is always a residual entropy at the static transition Trsb- 
As shown by Gross, Kanter and Sompolinsky [16] in the limit p — > oo the statics of the 
Potts model converges to the Random Energy Model (REM) [17, 27]. We will see that 
the convergence of the statics of the Potts glass model to the REM when p — > oo is 
very slow (logarithmic in p). Surprisingly, we will see that also the dynamics converges 
logarithmically with p to a fully frozen dynamics but even more slowly than does the statics. 
For all practical purposes, i.e. for reasonable values of p, the system is never fully frozen. 
In addition we will see that, for p > 4, the dynamical transition (also called the glass 
transition) is always very close to the static transition. This makes the glassy behavior of 
the Potts model very different from other models with a discontinuous transition in the 
order parameter like, for instance, the p-spin interaction Ising spin-glass model where the 
static and the dynamic Edwards- Anderson order parameter increase relatively fast with p. 

This partial freezing which occurs for the mean-field Potts glass has to be compared 
with deterministic models without quenched disorder and disordered spin-glasses. On the 
one hand, models like the low autocorrelation binary sequences [1] , fully frustrated lattices 
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in the mean-field limit [6] or discrete matrix models [5] display a stronger glass transition 
because qa ~ 1. In those cases there is no quenched disorder and frustration alone is the 
responsible to the existence of the glassy phase. The origin of the frustration is purely 
dynamical and self-induced by the dynamical process [3]. On the other hand, disordered 
models like the p-spin glass model [9] or the random orthogonal model (ROM) [2] do have 
a stronger freezing at the dynamical transition temperature. The presence of uncorrelated 
quenched disorder (i.e. JijJki = Jij Jm) in the Potts glass, and also in the SK model, has 
the effect of softening the discontinuous transition. In some sense frustration corresponds 
to the existence of some constraints on the different values of the quenched couplings J^. 
This is what happens in the random orthoghonal model where the form an orthogonal 
random matrix. Also in the case of the p-spin glass model the quenched disorder variables 
Jhi 2 ...i p tend to frustrate the system much more than in the Potts case. In the Potts glass 
case, for a fixed value of the number of states p, the number of quenched variables goes 
like iV 2 (N is the number of sites) while in the p-spin glass model this number increases 
much faster with iV for p > 2 (like JV!/(iV - p)\ ~ N p for finite p). 

The work is divided as follows. In section 2 we introduce the model and we write closed 
expressions for the free energy at first order of replica symmetry breaking. In section 3 
we solve numerically the static equations at one step of replica symmetry breaking and we 
determine the static and the dynamical transition. Section 4 compares the predictions to 
the numerical simulations. Finally we present the conclusions. 



2 Static replica equations for the Potts glass 

The Potts glass model is defined by the random Hamiltonian 

ft = -P^-^<W, > (!) 

i<j 

where p is the number of states and the variables a can take the values 0, 1, ..,p — 1. The 
sum is extended over all N( ^ N ~ 1 '> pairs in the lattice and iV is the number of sites. The 
couplings Jij are randomly distributed with mean ^ and variance In order to solve 
this random model we apply the replica method in order to compute the free energy / 

loff Z n 

Pf = lim n .o^ (2) 

where n is the number of replicas and the overline means average over the disorder. Per- 
forming the usual transformations (avergaging over the disorder, decoupling the sites and 
introducing auxiliary fields) and using the identity (a, b — 1, ..,n are replica indices) 
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p-1 

6 crfa b = J2 S KrS a b , (3) 

r=0 



one gets the following result 1 



Z"} = J dm r a dQ r a le- NA[m ' Q] , (4) 
where r, s — 0, — 1 denote the Potts states. The function A[m, Q] is given by 

A[m,Q] = np2{ \~ P) +^Uo + P -^) E«) 2 + ^EE(«) 2 -log Tr.e^M (5) 
with 



-(^0 + ^ ) E m a (P^„r- 

P Z a,r 

The stationary saddle point equations read 



P z a,r P a <b rs 



Q2 = ((p5* a r-l)(p5* b s-l)) (7) 

where the mean (...) is evaluated over the effective Hamiltonian in eq.(6). The order 
parameters m and Q satisfy the constraints 

E< = o 

r 

E« = . (8) 

r 

In the particular case p = 2 with Q r a f s = —Q a b, Q"b = Qab one recovers the solution 
for the SK model [8]. It can be shown that ferromagnetic order is always preferred for 
p > 2 for low enough temperatures. An estimate T# for the temperature Tp below which 
ferromagnetic order appears is given by the following condition [19] 

4+^ = 1 ■ (») 

In the special case J = the ferromagnetic transition appears below T = 1 for p < 4 
and above that temperature for p > 4. Our main interest in this work is the study of 
1 Alternatively one can use the simplex representation [18] 
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the spin-glass transition. In order not to observe the ferromagnetic transition it will be 
necessary to add an antiferromagnetic coupling in case p > 4. This will be discussed further 
at the end of this section. The spin-glass solution is given by m r a = 0. This means that all 
different p states are equally populated. The saddle point equations are independent of Jo 
and the replica symmetric solution in this case is given by 



Q2 = ?(p-i) ■ (io) 



Substituting this result in eq.(7) we obtain 

Pf='^r(l-p)(l-q) 2 - / II (^= e_ ") l0 §(E exp((3(qp)->y r )) . (11) 

The high-temperature result q = gives the free energy /, the internal energy u and 
the entropy s, 



„ = Mz& - log(p) 

« - ^ (12) 
3 2 (1 -p) 

s = \ +log(p). (13) 

Because the entropy has to be positive one finds that the replica symmetric solution 
breaks down, at least above or equalt to 



.41og(p), 

It has been shown [20] that there is a continuous phase transition at T c = 1 for p < 6 
which is unstable for p > 2. This transition ceases to exist above p = 6 and cannot be 
found within the replica symmetric hypothesis. 

It is necessary to break the replica symmetry. By expanding the free energy eq. (5) close 
to T c = 1, Gross, Kanter and Sompolinsky [16] have found two different regimes according 
to the value of p. In both cases the correct solution is given by one step of breaking. In the 
region 2.8 < p < 4 the transition is continuous. The breaking parameter m is E ^ at the 
transition temperature T c — 1. At low enough temperatures the entropy of the one-step 
solution becomes negative and a continuous breaking is then necessary. In the regime p > 4 
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the transition is discontinuous in Q and the breakpoint parameter m is equal to 1 at the 
transition temperature T c > 1. Cwilich and Kirkpatrick [21] have shown that this one step 
solution is always stable for p > p* = 2.82 below but close to T c . 

At first order of replica symmetry breaking we subdivide the n replicas into ^ blocks. 
Each block contains m replicas [22]. The order parameter Q™ b takes a certain value when 
both replicas a, b belong to the same subblock and it is zero when both indices belong to 
two different subblocks. More explicitely, if K denotes a subblock, we impose 



Qlt " = -Q K ie 4 Qat " = (otherwise) 
Q r ; b = -q(p-l) (a,beK), QZ = (otherwise) (15) 



We obtain the result, 



B 2 B 2 B 2 B 2 

Pf=<j(l-p) + ^(m-l)(p-l)q 2 + ^q(p-l) + P -qm- 



1 roo P A v y 2 P 

~ log / II (-W= e "" f ) (E exp(P(qp)-iy r )) m (16) 
The corresponding saddle point equations are 

dq dm ^ ^ ^ 

which determine the correct solution. It is possible to solve pertubatively these equations 
in three different cases: 

• Expanding around p = 4 since in this limit case the transition is quasi-continuous 
[12, 13]. This technique has be applied also in case of the p-spin model [11]. 

• Expanding around T = 1 using the so-called effective approximation for the free 
energy eq.(5) up to order Q A [21]. 

• Solving the limit p — > oo. In this limiting the model converges to the random energy 
model [23]. For recent work see [24] 

We are interested in the glassy behavior of the Potts model. Our approach will be to 
numerically solve the equation (17). This is the purpose of the next section. 

Some comments are in order regarding the existence of the ferromagnetic transition. 
We said previously that the system orders ferromagnetically at low enough temperatures. 



7 



The temperature Tp below which the system orders ferromagnetically is smaller than Tp 
with Te < 1 for p < 4. Also for p < 4 the spin-glass transition appears at Trsb = 1- 
This means that in the regime p < 4 the spin-glass transition occurs at a temperature 
Trsb larger than the temperature Tp at which ferromagnetic order sets in. On the other 
hand, for p > 4 the spin-glass transition Trsb occurs at a temperature greater than 1 but 
smaller than Tp. In order that Tp < Trsb it is necessary to introduce a negative value 
for J . In our numerical simulations we have chosen J = -^p in case p > 4 and Jq = 
for p < 4. In this way the spin-glass transition occurs at a larger temperature than the 
ferromagnetic ordering. Now, let us suppose that we perform a dynamical process of the 
system in which the temperature is slowly decreased starting from the high-temperature 
phase. We think that, once the system has entered the metastable glassy phase, then it 
remains trapped in this phase for a time which diverges exponentially with the size of 
the system. Consequently, the system is unable to see the ferromagnetic transition which 
occurs at a lower temperature. Two reasons reinforce this observation: 

• Numerical studies of the case Jq = for p > 4 (see section 4) show that the ferro- 
magnetic transition occurs at a temperature Tp smaller than Tp (eq.(9)). Then, for 
a generic negative value of Jo we can expect the ferromagnetic transition to appear 
at a temperature much lower than Trsb- 

• We can expect there exists a glass transition associated to the static spin-glass one 
which probably occurs at a temperature Tq larger than Trsb- 

In summary, by choosing J as indicated above, the ferromagnetic transition tempera- 
ture will always be smaller than the freezing temperature at which the spin-glass ordering 
appears. Only the case p = 4 could be a little tricky because the estimate for the ferro- 
magnetic transition Tp and the spin-glass transition coincide, but even in this case we have 
not observed in the numerical simulations a strong magnetic ordering. 

We should note that there are very few works devoted to the study of the ferromagnetic 
behavior in the mean-field Potts glass and we think it would be interesting to investigate 
it. 

3 The static and the dynamical transition 

In this section we are going to solve numerically eq.(16). As is usual in spin-glass theory 
we have to maximize the free energy as a function of q and m. We face the problem of 
computing the p-dimensional integral 
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POO P An, „2 P , 

1 = / II (^= e ^) (E exp(/5(gp)^ r )) m (18) 

Because the solution of the replica equations involve a maximization in the (q, m) plane 
it is mandatory to compute I with relatively high precision. Because J is a p-dimensional 
integral it can be computed with the usual techniques only for p not too large. The most 
easy thing one can do is to divide the p-dimensional space into small squares and use the 
Simpson algorithm or a similar one. The computation time grows as a power of p. This 
makes the calculation practically impossible for p > 3. However, if one exploits the fact that 
the integration argument is invariant under the permutation of the p indices r = 0, ..,p — 1 
then the integration region can be reduced to the hyperplane y\ < y<i < ... < y p where 
the y r denote the coordinates of one point. In this way one can gain a factor p\ in the 
computation time and we have been able to solve with good enough precision up to p = 7. 
These preliminary methods should be considered as checks for the main computation. 

We have been able to reduce the p-dimensional integral to a two dimensional integral. 
In this way the problem is completely resoluble, at least numerically. We use the identity, 

j f-OO 



[ dxx~ m e~ Ax (19) 
Jo 



r(i -m) 

where T(x) is the Gamma function that is well defined for x > 0, i.e. m < 1 as is the case 
once the analytic continuation n — > (n is the number of replicas) has been done. 

We decompose the integrand in eq.(18) as a product of two terms A * A m ~ x with A 
given by, 

p i 

A = Yl exp(P(qp)*y r ) (20) 

r=l 

Applying eq.(19) and using the fact that the integrand in eq.(18) is invariant under 
permutation of the indices we get the final result, 



(§_qp_\ 

pe y 2 > 



r(i -m) 

where w(x) is given by 



[°° dxx- m (w(x)) p - 1 w(xe^ E ) (21) 
Jo 



/OO 1 
— ==exp(-— - x exp(p(qp)iy)) (22) 
-oo \/Z7T 4 



dy y 2 

-oo y/2n 

The integral over x is well defined and free of divergences. However one has to be 
careful evaluating the integrand close to x — 0. We have been able to maximize the free 
energy and completely solve the static replica equations up to p = 40. 



9 



The results are shown in figures 1 and 2 where we plot the variational parameters q 
and m as a function of T. We plot the solutions for the cases, p — 3, 5, 7, 10, 20, 40. The 
transition temperature also grows with p. The solution of the integral eq.(21) presents some 
problems of precision at very low temperatures and also close to the transition temperature 
where it is difficult to precisely determine the value of the discontinuity. A more precise 
way to compute the critical temperature and the discontinuous jump of q will be presented 
below. It is interesting to note how much slow is the convergence to the random energy 
model as p increases. When p increases the value of q at the transition point grows very 
slowly. In fact, in the limit p — > oo, the value of q converges to 1 and the entropy is zero at 
the transition temperature. Using eq.(14) we obtain that the critical temperature grows 
like T of eq.(14). This result was already noted in [16]. 

We have already observed that at very low temperatures the entropy of the one-step 
solution becomes negative (of order 10~ 2 ). Continuous breaking is necessary (as noted in 
[16]) but we have not studied this type of solution. It is not clear to us if any effect of this 
new transition could be observable in a numerical simulation. 

We want to show now a more precise computation of the critical temperature T r $b and 
the glass transition Tq. From the dynamical point of view an instability in the dynamical 
equations appears at a temperature Tq above the static transition T^sb- Using the static 
approach, this dynamical temperature T G can be determined computing the smallest eigen- 
value in the stability matrix. The vanishing at Tq of this eigenvalue, sometimes called the 
replicon, corresponds to the marginality condition [25]. In principle, this condition cor- 
rectly determines the dynamical or glass transition. Anyway it is not clear if it is the 
correct description of the dynamical behavior in the low temperature phase. This condi- 
tion has been numerically solved in the random orthoghonal model and it has been shown 
that it correctly describes the dynamical energy below the glass transition and not too low 
temperatures [2]. It also describes correctly the glass transition in case of deterministic 
models. The interested reader is referred to [2] for more details. In order to determine the 
glass transition for the Potts case we should compute the stability matrix of the problem. 
This is an involved task (which has beeen done by Cwilich and Kirkpatrick close to T c [21]) 
and we will follow a different strategy (already noted by Cwilich and Kirkpatrick but not 
fully explained). It can be shown that in the limit m — > 1 the replicon eigenvalue coincides 
with the longitudinal eigenvalue. This result can be shown using the exact expressions 
for the spectrum of the stability matrix which have been reported in the literature at first 
order of replica symmetry breaking [26] . From the stability analysis results of Cwilich and 
Kirkpatrick this can also be directly tested in the Potts glass case. Consequently, in order 
to determine the dynamical transition, suffices to impose the marginality condition for the 
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longitudinal fluctuations. 

We expand the free energy eq.(16) around m — 1, 

Pf = \P 2 (1 ~ P) ~ log(p) + (m - 1) Q/3 2 (p - l)q 2 + \(3 2 q{p + 1) + log(p) - J 2 ) (23) 
where the integral 7 2 is given by, 

I 2 = exp(-^) / J] (-#=e-f ) l og (^ exp((3( qP ^y r )) (24) 

^ J-oo r=1 y27T r=1 

For m = l eq.(23) reduces to the high-temperature free energy which is independent of 
q. More generally, we can write the free energy as 

f = f 0+ (m-l)f 1 + O((m-l) 2 ) (25) 

where fo is independent of q. This general expansion locates the static and the dynamic 
transition. For the static transition we look for a temperature at which there is a solution 
qRSB which satisfies 



d£\ _ (dfi 



Ui) q = qRSB = (26) 

For the dynamical transition the stability is marginal and the second derivative of / 
respect to q vanishes, 



dq 2 )q=q G {dq 2 )q=q G ' ^ 

The last equations correspond to the case in which an extremal solution of eq.(23) with 
qc 7^ dissappears. It is then clear that the dynamical transition temperature is always 
higher than the static one. We have solved these equations for different values of p. Now 
we face the problem of computing the p-dimensional integral It can be reduced to a 
two dimensional integral using the representation, 

log(l + A) = / — e~ x (l - e~ Ax ) . (28) 
Jo x 
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and taking 



A = (E exp(/3(qp)*y r )) - 1 (29) 

r=l 

we obtain the result 

/ 2 = r — e - x {l-e x w(xe^)w p - 1 (x)} . (30) 

Jo x 

with the same function w as given in eq.(22). We have solved equations eq.(26) and 
eq.(27) for different values of p. Our results are summarized in Table 1. We find for each 
value of p two temperatures. One is Trsb and corresponds to the static transition with 
the discontinuous jump qRSB- The other one is T G and corresponds to the dynamical 
transition with the discontinuous jump q G . Our results for the static transition are in 
agreement with those found with the previous analysis using the maximization procedure 
for the free energy This is a check of our procedures. Moreover this analysis provides a 
much more precise determination of the values of qnsB, qG and the transition temperatures. 

The results we have found are also in agreement with those reported by Cwilich and 
Kirkpatrick, the only difference is that all their computations are perturbative whereas 
ours are exact. As was obtained in [12] and [21] one finds that q G /qnsB = 3/4 for p close 
to 4. Looking at table 1 the reader can observe that the ratio q G /qRSB stays so close to 
3/4, even for large values of p, that one is tempted to conclude that this is exact at all 
orders in perturbation theory. Our numerical precision to solve the equations (26), (27) is 
good enough to exclude this possibility. From the results shown in the Table 1 it is clear 
that the convergence to the p — > oo limit is very slow. Fortunately, our numerical program 
which solves the equations (26) and (27) is enough accurate to explicitly show this slow 
convergence even for exponentially large values of p. We have solved the full equations up 
to p — 10 6 . The results for qRSB and q G as a function of j^^y are shown in figure 3. 

Furthermore in the Potts case the ratio T g /T RS b grows very slowly with p being always 
smaller than 1.13 up to p = 10 6 . The proximity of the temperatures makes it difficult to 
discern one from the other in numerical simulations. This proximity of the static and 
the dynamic transition temperatures is very probably related to the small value of the 
dynamical order parameter q G for large values of p. From these results we expect the 
glassy behavior of the Potts glass to be very different from other disordered spin-glass 
models. 

For instance, in case of the p-spin interaction spin-glass model we have also solved the 
equations corresponding to eq.(26) and eq.(27). We have found that both the static order 
parameter q^sB and the dynamical q G converge to 1 in the limit p — > oo much faster than 
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the Potts case in agreement with theoretical expansions around the p — > oo limit [24] . For 
p = 3 (the smallest value of p compatible with a discontinuous transition) one finds in the 
p-spin model, 

q RSB ~ 0.81 (Trsb ^ 0.65) 
q G ~ 0.68 (T G ^ 0.68) . (31) 

For this particular model, the ratio qa/qRSB tends to 1 in the limit p — > oo and the 
ratio T G /T RSB increases with p much faster than the Potts model does (for p = 10 we find 
qc/qRSB ^ -97 and T G /T RSB ~ 1.38) 

In the next section we shall compare all these predictions with Monte Carlo numerical 
simulations. We will see that the Potts glass transition is always present but it is far 
from being a complete thermodynamic freezing as happens in models where frustration is 
stronger (see, for instance the random orthogonal model [2]). Before showing our Monte 
Carlo results for the spin-glass transition it will be interesting to present some results on 
the ferromagnetic ordering that takes place in the Potts glass. The problem of the existence 
also of a real glass transition above the ferromagnetic transition still remains open. Our 
main interest is to show that if one does not introduce an antiferromagnetic coupling then 
the ferromagnetic ordering takes place even though the transition temperature is well below 
that given by eq.(9). 

4 Monte Carlo tests of the glass transition 

In order to simulate the Potts glass we have considered the Hamiltonian, 

h = -Y / Mp^ j -i) (32) 

i<j 

The Jij are distributed with mean Jq/N and variance equal to For computational 
reasons we have chosen a binary distribution where the can only take the values ±77^- 
The only difference between the Hamiltonians of eq.(32) and eq.(l) is a constant which 
vanishes in the thermodynamic limit. We have chosen this second version because we have 
found that the addition of the constant strongly reduces the sample to sample fluctuations 
in the high-T region. This should not make too much difference for small values of p but 
is crucial for large values of p. All simulations implement the Metropolis algorithm with 
random updating. 
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The results we present in the next subsections correspond to annealings in which we 
compute the main thermodynamic observables. Starting from the high-temperature region 
the temperature is progressively decreased. Statistics is collected at each temperature and 
the time we stay at each temperature is the same for all temperatures during the cooling 
procedure. We have computed the internal energy, the magnetization of the different in- 
states and the associated dissipative quantities, i.e. the specific heat and the p different 
magnetic susceptibilities corresponding to each one of the p-states. The specific heat and 
the magnetic susceptibility of one of the p states is computed measuring the fluctuations of 
the internal energy and the magnetization (see eq.(7)) respectively. Typically we performed 
several thousands of Monte Carlo sweeps at each temperature. We have to call the attention 
of the reader that our results are dependent on the time schedule of the annealing only 
for very large values of p (i.e. where the finite-size corrections are large). Otherwise, one 
cannot observe a sensible dependence of the different quantities on the time the system 
stays at each temperature during the cooling procedure. At least, this dependence is of 
the same order as that arising from the sample to sample fluctuation. Obviously we have 
performed annealing schedules as large as possible within our computing capabilities. We 
will eventually show the dependence on the annealing time in the large p case. 

4.1 Ferromagnetic ordering with Jo = 

When J = the system orders ferromagnetically. We have investigated the ferromagnetic 
ordering for p = 10. This value of p is in the regime (p > 4) where the ferromagnetic 
transition is expected to appear at a temperature higher than the spin-glass transition. 
From eq.(9) we expect ferromagnetic order to be present below T E = 4. Figure 4 shows the 
internal energy as a function of the temperature compared to the energy of the spin-glass 
phase and the high-temperature result eq.(13). The energy is lower than that corresponding 
to the spin-glass solution. Figure 5 shows the result for the magnetic susceptibility averaged 
over the different p = 10 states. The first well defined peak is at Tp = 2.3 which corresponds 
to the temperature at which the energy of the system departs from the high-temperature 
result (dashed line in figure 4). The specific heat also shows a peak at that temperature. 
The reader will immediately recognize that temperature as the ferromagnetic ordering 
temperature Tp which is much lower than the estimate Tp. It is interesting to note in figure 
5 the emergence of further peaks at lower temperatures. These are a sign of new transitions 
in the ferromagnetic phase. From the measurements of the magnetization we have observed 
that at the transition temperature there is one state which acquires a magnetization greater 
than zero which means that one state is macroscopically populated. We also have observed 
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that the other peaks at lower temperatures correspond to the emergence of new states which 
start to be macroscopically populated. 

We have also studied the zero temperature ground state following a steepest descent 
procedure. We have searched for stable configurations against one-spin flip (for a general- 
ization to stability against n-spin flips see [28]). Starting from a random initial configura- 
tion we sequentially move on the lattice selecting (among the p — 1 possibilities) the state 
of the variable a(i) which releases the largest amount of energy. In this way the system 
reaches a metastable state that should be magnetized if the ground state is ferromagnetic. 
We have repeated this procedure several times saving the energy and the magnetization of 
the final configurations. Figures 6 and 7 show the distribution probability of the energies 
and the magnetization eq.(7) of the stable configurations against one-spin flip movings for 
the same model J = 0, p = 10 with N = 100. Figure 6 shows that the energies of this 
class of metastable states are distributed very similarly to the form predicted for the SK 
model [29]. This is a consequence of the glassy nature of the ferromagnetic phase. We 
have verified that the minimum energy found by the algorithm is higher than the energy 
obtained doing a slow cooling starting from the high-temperature phase. This is a proof 
of the glassy nature of this phase. The ferromagnetic nature (but glassy) of the phase is 
explicitely shown in figure figure 7. The value of m ranges from m = — 1 to m =p - 1 = 9. 
The peak at m = — 1 is consequence of the fact that only some of the p — 10 states are 
populated. 

4.2 The continuous transition (p = 3) 

The case p = 3 is indeed very similar to the SK model (p — 2). Because the transition 
is continuous the system relaxes very close to the true energy. As mentioned in section 
2, it now suffices to take J = 0. In this way the ferromagnetic transition lies well below 
the spin-glass transition. In fact, we have not observed any tendency of the system to 
be magnetized at low temperatures. Figure 8 shows the internal energy as a function of 
the temperature along with the one-step solution and the high-temperature result eq.(13). 
Below the critical temperature Trsb = 1 the data departs from the prediction. Precisely 
at T = 1 the specific heat and the magnetic susceptibility have a cusp. Similar results for 
the internal energy were obtained for p — 4. 

4.3 The discontinuous transition 

To investigate the discontinuous spin-glass transition we have chosen J = for p > 
4. In this way the system first enters the metastable glassy phase in which there is no 
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ferromagnetic ordering. In all our simulations we have observed that this is what happens 
and that there is no tendency for the ferromagnetic domains to grow as the temperature 
is decreased. At high temperatures the size of the domains (i.e. the fraction of sites of the 
lattice which are in the same state) is 1/p. This is true down to very low temperatures 
where in the worst case the size of the domains increase approximately ten percent. To 
make any tendency to the ferromagnetic ordering completely dissapear we can increase 
the intensity of the antiferromagnetic coupling. This is only possible if p is not too large 
because otherwise finite-size corrections (an consequently finite-time effects) considerably 
increase. In the regime of large values of p one can neglect finite-size effects only if p <C N. 
This is a problem of the simulations in the large p regime (the region where the glass 
transition can be clearly appreciated). We will explicitly show the time dependence of the 
results of annealing for p = 20. We divide our results for the discontinuous transition in 
two parts, the small and the large p > 4 regime (corresponding to the magnitude of the 
finite-size corrections). 

4.3.1 The small p > 4 regime 

We have measured the internal energy as a function of the temperature for cases p = 5, 10. 
In this regime we have observed that the results do not vary too much depending with 
the time schedule of the cooling procedure. Comparison with theory is shown in figures 
9 and 10. For these small values of p the dynamical transition practically coincides with 
the static one. Comparing to the previous continuous case p = 3 we see that the energy 
in the low T region is very close to the predicted one for p = 5 and remains slightly 
above the expected one in case p = 10. This is the glassy phase where the system remains 
trapped making excursions between several metastable states of similar energy but without 
reaching the static phase of slightly lower free energy. The difference in free energy (and 
energy) between the static phase and the metastable glassy phase is small for p = 5 and 
increases with p. It is important to note that the energy we are computing is purely 
dynamical. For p < 4 this difference of free energy does not exist. This does not mean 
that the system relaxes to the true ground state energy in an annealing process (see figure 
8). In fact, for a continuous transition we expect the system should relax to the static free 
energy at a finite temperature very slowly (as a power law) very similarly to the relaxation 
of the remanent energy or the remanent magnetisation in the SK model [30, 31]. For a 
discontinuous transition the relaxation of the free energy takes place also very slowly but 
to a dynamical value higher than that predicted by the static approach. We have also 
computed the specific heat and the magnetic susceptibility. They display a cusp located 
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approximately at the static transition (and, because of its proximity, also the dynamical 
transition). The reader may be a little puzzled by the data shown in figure 10 because the 
dynamical energy departs form the high-T behavior at a temperature higher than the glass 
transition temperature. We think this occurs because the static and the glass transition 
temperatures are very close to each other. We will return to this point in the conclusion. 

4.3.2 The large p > 4 regime 

Finite-size corrections are important and one has to simulate large sizes in order to reduce 
these effects. We present the results of annealings for p = 20,40 in figures 11 and 12. We 
decided to simulate the Gaussian model instead of the binary ± J one in order to reduce 
finite-size corrections. As p increases the finite-time effects also increase and we have found 
a clear dependence of our results on the cooling procedure. Figure 11 shows simulation 
results for p = 20 for N = 2000 and two different cooling procedures. The simulation 
results show a drift with the annealing time. For p = 40 (figure 12) we show simulations of 
two different sample realizations. Since sample to sample fluctuations increase with p, the 
relative magnitude of the fluctuations of figure 12 should be considered as an upper bound 
for the previous figures with smaller values of p. Also finite-time effects are large for p = 40 
. From figures 11 and 12 it can be appreciated that also in this case the energy departs 
from the high-T result at a temperature higher than predicted for the glass transition. 

Glassy effects are much more pronounced in the large p-regime, the dynamical energy 
being larger than the static one. All dissipative quantities show a cusp very close to the 
dynamical transition. Even though the static and the glass transition are very close one to 
the other the fact that the energy of the sytem is much higher than the static one (when 
approaching the glass transition) is a proof that the system has entered the glassy phase. 

5 Conclusions 

We have studied the glassy behavior of the mean-field Potts glass. We have been able to 
numerically solve the static equations at first order of replica symmetry breaking. Also 
we have introduced a simple method, already observed by Cwilich and Kirkpatrick [21], 
which allows a full computation of the static and the dynamical or glass transition and the 
associated Edwards-Anderson parameter. 

We have numerically computed the parameters of the transition for different values of 
p. We observe that the Edwards- Anderson parameter at the glass transition qc increases 
logarithmically with p. The ratio of the static and glass temperature is smaller than 1.13 
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up to p = 10 6 . The situation is very different from other disordered models such as the 
p-spin Ising model. For that model the dynamic transition temperature is much higher 
than the static one. 

All our numerical results seem to indicate that the dynamical transition takes place at 
a temperature higher than that predicted by the theory. But this is due to the proximity 
of the dynamical transition to the static one. If the dynamical transition temperature were 
much larger than the static one then we would expect that the energy departs from the 
high-T result precisely at the dynamical temperature. This is indeed the situation one 
observes in low autocorrelation models [1, 3, 7] , in the random orthogonal model [2, 6] 
and discrete matrix models [5] . We would also expect this situation for the p-spin Ising 
model in the regime of not too small p, but unfortunately we are only able to perform 
simulations for p = 3. 

For large values of p (like p = 20,40) we have observed a clear dependence on the 
time spent during the cooling procedure . The origin of the finite-time effects is related 
to the finite-size effects we also observe in this regime. We expect that simulations for 
larger sizes should give results nearly independent of the annealing time leaving only a 
small thermalization time effect close to the glass transition where critical effects begin to 
be important. We interpret this effect in the following scenario. 

There are two characteristic relaxation times in the system. The first time tq diverges 
as the dynamical transition is approached, the other one r s diverges as the static transition 
is approached. Because T G is so close to Trsb the systems feels the static low temperature 
phase very close to the dynamical transition temperature. Above the glass temperature 
we have t g ~ t s which is certainly large if the system is entering the low temperature 
phase. Because the characteristic time scale tq increases very fast only very close to 
the dynamical transition temperature then we expect that close to Tq the correlation 
time r s will set the characteristic time scale above which our simulation results should 
be time independent. Only for times larger than r s (which we are not able to reach in 
our simulations) the system would behave as dynamics predicts. It is then clear that all 
our simulation results are smeared by the static relaxation time r s . In the other models 
mentioned in the previous paragraph the dynamical transition temperature is much larger 
than the static one. Approaching the dynamical transition the system is in the high 
temperature phase where the relaxation time is very small. Consequently, tq grows very 
much only very close to T G (very probably diverges like r G ~ (T — T G )~ 7 where 7 = 2 the 
typical value for mean- field models [7]) and the system departs from the high-T result very 
close to that temperature. 

Because the Potts model is only partially frozen at the glass transition this is a model 
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appropiate for study of the dynamics in the metastable glassy phase. We expect that 
jumps among states could be seen numerically without special effort. For models without 
disorder, the system freezes quickly at the glass temperature and more involved numerical 
techniques are needed in order to allow the system to change state [32]. 

We would also like to draw attention to the interest of studying the ferromagnetic 
ordering for zero mean coupling. It would be interesting to understand the static as well 
as the dynamical behavior in that case. 
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Figure caption 



Fig. 1 The one-step breaking parameter q as a function of the temperature. From left to 
right: p = 3, 5, 7, 10, 20, 40 

Fig. 2 The one-step breaking parameter m as a function of the temperature. The different 
lines that intersect the upper horizontal axis m = 1 correspond from left to right to: 
p = 3,5,7,10,20,40 

Fig. 3 The static and the dynamic Edwards- Anderson parameter as a function of j^^y- 
They increase logarithmically with p. The dots are for the static value, the crosses 
for the dynamical one. 

Fig. 4 Energy versus temperature for the case p = 10 with Jo = 0. The continuous line 
corresponds to the one step spin-glass solution and the dashed line is the high-T 
result. The ferromagnetic transition is close to 2.3. Simulation results are for one 
sample and N = 1000. 

Fig. 5 Magnetic susceptibility versus temperature for the case p — 10 with J = 0. The first 
peak appears at T F ~ 2.3. 

Fig. 6 Probability distribution of the energy of the ground states for case p = 10 and J = 
and N = 100. 

Fig. 7 Probability distribution of the magnetization associated to the zero-temperature 
metastable states. The singularity at m = —1 means that only some states are 
populated. The same parameters as in figure 6. 

Fig. 8 Energy versus temperature in case p = 3, J = 0. The continuous line is the one- 
step solution. The dashed line is the high-T result. The transition temperature is 
T RSB =T G = 1. The full dots are for one sample and N = 2000 

Fig. 9 Energy versus temperature in case p = 5, Jq = —1. The continuous line is the one- 
step solution. The dashed line is the high-T result. The transition temperatures are 
T RSB ^ 1.0091 , Tq ~ 1.01. The full dots are simulation results for N = 1000. 

Fig. 10 Energy versus temperature in case p = 10, Jo = —3. The continuous line is the 
one-step solution. The dashed line is the high-T result. The transition temperatures 
are T RS b — 113 ,T G ~ 1.14. The full dots are simulation results for iV = 1000. 
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Fig. 11 Energy versus temperature in case p = 20, J = —8. The continuous line is the one- 
step solution. The dashed line is the high-T result. The transition temperatures are 
Trsb — 1.36,Tg ~ 1.39. The crosses and the full dots correspond to the Gaussian 
model with N = 2000 and two cooling procedures (the simulations with crosses are 
10 times larger in simulation time than the dots). 

Fig. 12 Energy versus temperature in case p = 40, J = —18. The continuous line is the 
one-step solution. The dashed line is the high-T result. The transition temperatures 
are T RSB ^ 1.71,T G ~ 1.76. The crosses and the full squares correspond to the 
Gaussian model with N = 1000 and two different samples. 
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